STABILITY OF THE STOCHASTIC HEAT EQUATION IN ^([0,1]) 

NICOLAS FOURNIER AND JACQUES PRINTEMS 

Abstract. We consider the white-noise driven stochastic heat equation on [0, oo) X [0, 1] with Lipschitz- 
continuous drift and diffusion coefficients b and a. We derive an inequality for the L 1 ([0, 1])-norm of 

>»»2 the difference between two solutions. Using some martingale arguments, we show that this inequality 

provides some a priori estimates on solutions. This allows us to prove the strong existence and (partial) 
uniqueness of weak solutions when the initial condition belongs only to L 1 ([0, 1]), and the stability of 
the solution with respect to this initial condition. We also obtain, under some conditions, some results 
concerning the large time behavior of solutions: uniqueness of the possible invariant distribution and 

VO , asymptotic confluence of solutions. 
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1.1. The equation. Consider the stochastic heat equation with Neumann boundary conditions: 

{d t u(t,x) = d xx u(t,x) + b(u(t,x))+a(u(t,x))W(t,x), i > 0, x G [0, 1], 
u(0,a:) = uo(a:), see [0,1], 
d x u(t,0) = d x u(t,l) = 0, t>0. 

Here b, a : R H> R are the drift and diffusion coefficients and uq : [0, 1] i— > R is the initial condition. We 
write formally W(dt, dx) = W(t, x)dtdx, for W(dt, dx) a white noise on [0, oo) x [0, 1] based on dtdx, see 
Walsh [13] . We will always assume in this paper that b, a are Lipschitz-continuous, that is for some C, 
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1. Introduction and results 



Our goals in this paper are the following: 

• prove a strong existence and (partial) uniqueness result when the initial condition uq only belongs to 
L 1 ([0, 1]) and some stability results of the solution with respect to such an initial condition; 

• study the uniqueness of invariant measures and the asymptotic confluence of solutions. 

We will investigate these two points by using some a priori estimates on the difference between two 
solutions u, v, obtained as a martingale dissipation of the L ([0, l])-norm of u(t) — v(t). 

Let us mention that our results extend without difficulty to the case of Dirichlet boundary conditions 
and to the case of the unbounded domain R (with uq G L 1 (R)). 

This equation has been much investigated, in particular since the work of Walsh (T3] . In [TJ] , one can find 
definitions of weak solutions, existence and uniqueness results, as well as proofs that solutions are Holder- 
continuous, enjoy a Markov property, etc. Let us mention for example the works of Bally-Gyongy-Pardoux 
[1] (existence of solutions when the drift is only measurable), Gatarek-Goldys [7] (existence of solutions 
in law), Donati-Pardoux (comparison results and reflection problems), Bally-Pardoux (smoothness of the 
law of the solution), Bally-Millet-Sanz [3] (support theorem), etc. Sowers [12], Mueller [9] and Cerrai [4] 
have obtained some results on the invariant distributions and convergence to equilibrium. 
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1.2. Weak solutions. We will consider two types of weak solutions, which we now precisely define, 
following the ideas of Walsh [T3j . When we refer to predictability, this is with respect to the filtration 
(-7 r t)t>o generated by W, that is T t = a(W(A), A G B([0,t] X [0, 1])). 

We denote by L' p ([0, 1]) the set of all measurable functions / : [0, 1] h-> R such that ||/||lp([o,i]) = 
(£ \f(x)\Pdx)^ < 00. 

Finally, we denote by Gt(x, y) the Green kernel associated with the heat equation dtu = d xx u on R + x [0, 1] 
with Neumann boundary conditions, whose explicit form can be found in Walsh [13} . Here we will only 
use that for some Ct, for all i,t/£ [0, 1], all t e [0,T], see [T3] . 

(2) 0<G t (x,y)<^e-^ 2 / 4t . 

Definition 1. Assume {%), and consider a R-valued predictable process u = (u(t, aO)t>o,xe[o,i]" 
(i) For uq G L ([0, 1]), u is said to be a weak solution to {7J) starting from uq if a.s., 

(3) forallT>0, sup ||w(*)||ii([o,i]) + / ll cr ( w (*))llz, 2 (fo i]) dt < °° 

[o,T] Jo 

and if for all (p G C\ ([0, 1]) such that tp'(Q) — </?'(l) = 0, for all t > 0, a.s., 

(4) / u(t,x)ip(x)dx = / uo(a;)v ; '(a ; )rfa ; + / / o-(u(s,x))ip(x)W(ds, dx) 
Jo Jo Jo Jo 

ft /-l 

/ [u(s, x)ip"(x) + b(u(s, x))ip(x)]dxds. 
'o Jo 

(ii) For Uq bounded-measurable, u is said to be a mild solution to {!]] starting from uq if a.s., 

(5) for all T > 0, sup \u(t,x)\ <oo 

[0,T]x[0,l] 



(in 



d if for all t > 0, all x £ [0, 1], a.s., 



(6) u(t,x)= / G t (x,y)u (y)dy+ / G t _ s (x,y)[o-(u(s, S/))W(ds, dy) + b(u(s,y))dydsj. 

Jo Jo Jo 

Let us make a few comments. Recall that for (H(s, y)) s >o j/e[0 l] a R-valued predictable process, the 
stochastic integral L L H(s, y)W(ds, dy) is well-defined if and only if L L H 2 (s, y)dyds < oo a.s. 

• Thus ([3]) implies that all the terms in ((4]) are well-defined. Clearly, condition ([3]) is not far from minimal. 

• Next, ([S]) and (O imply that all the terms in (JSJ) are well-defined, but here |(SJ) is clearly far from optimal. 

When uo only belongs to i 1 ([0, 1]), we will only be able to prove that ([3]) holds. 

Let us finally recall that Walsh [T3] proved, under (H), that for any bounded-measurable initial condition 
Uo, there exists a unique mild solution u to ([1]), which is also a weak solution and which furthermore 
satisfies, for all p > 1, all T > 0, E[sup[ .t]x[o,i] l M (*' :z; )l P ] < °°- 

1.3. Existence and stability in L 1 ([0,1]). Our first goal is to extend the existence theory to more 
general initial conditions. 

Theorem 2. Assume {H). 

(i) For uo G L ([0, 1]), there exists a weak solution u to (QJ) starting from ug. 

(ii) This solution is unique in the following sense: for any sequence of bounded-measurable functions 

Uq : [0, 1] i — >• IR. such that lim n ||u l — uo||ii([o ] ii) =0, the sequence supr T i \\u n (t) — u(t)|| i iQo,ii) tends to 

in probability for any T. Here u n is the unique mild solution to Uty starting from Uq. 
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(Hi) For Uq,v G i 1 ([0, 1]), consider the two weak solutions u and v to Q\) starting from u and Vq built 
in (i). For all 7 G (0, 1), all T > 0, we have 



E 



T \ 7/2 ' 



sup |K*) - uO)H2i([o,i]) + I / lk(«(*)) - c(K*))IL=([o,i])* 



< C&,7,t||U0-W ||2i([0,i]), 



where Ch t ~/,T depends only on b, 7, T. 

(iv) Assume now that b is non-increasing. For uo, vq G L ([0, 1]), tei u, w fee t/ie two weafc solutions to (QJ) 

starting from uq and vq built in (i). For all 7 G (0, 1), we ftawe 



sup ||«(t) - v(t)\\l H[0 !])+(/ ||6(«(t)) - 6(«(«))||ii([o,i])* 

0,00) ' VJo j 

' H cr (wW)-CT(w(i))||| 2([0!l]) dtJ < C 7 ||uo-^o|lli([o,l]): 



where C 7 depends only on 7. 

Observe that this result contains a regularization property. For example if o~(z) = z, even if uo does not 
belong to L 2 ([0,1]), the weak solution satisfies ((3|) and in particular a(u(t)) = u(t) € L 2 ([0,1]) for a.e. 
t > 0. For the same reasons, the stability result (iii) provides a better estimate for a.e. t > than for 
t = 0. 

To our knowledge, Theorem [2] is the first result concerning L 1 ([0, 1]) initial conditions. Many works 
concern bounded-measurable (or continuous) initial conditions, see Walsh |13j . Bally-Gyongy-Pardoux 
PP, Cerrai [I]. Another abundant literature deals with the Hilbert case (initial conditions in L 2 ([0, 1])), 
see Pardoux [TU], Da Prato-Zabczyk [5], Gatarek-Goldys [?]• 

The present well-posedness result is quite satisfying, since the requirement that wo G -^ 1 ([0, 1]) is very 
weak and seems necessary for (j4]) to make sense. 

1.4. Large time behavior. We now wish to study the uniqueness of invariant measures. 

Definition 3. A probability measure Q on L 1 ([0,1]) is said to be an invariant distribution for (QJ) if 
for uq a ^ 1 ([0, l])-valued random variable with law Q independent of W, for u the weak solution to J7]) 
starting from Uq built in Theorem^ C(u(t)) — Q for all t > 0. 

We have the following result. 

Theorem 4. Assume (H), that b is non-increasing and that (<r, b) : R n- K 2 is injective. Then (QJ) admits 
at most one invariant distribution. 

To prove the asymptotic confluence of solutions, we need to strengthen the injectivity assumption. 



There is a strictly increasing convex function p : K + h-> M + with p(0) — such that 
forallr,zGR, \b(r) - b{z)\ + \a(r) - a(z)\ 2 > p(\r - z\). 



Theorem 5. Assume (H), that b is non-increasing and (I). 

(i) The following asymptotic confluence property holds: for uq,vq G i 1 ([0, 1]), for u, v the weak solutions 

to (Op starting from uq and vq built in Theorem @ 

a.s., lim \\u(t) - f (t)| |j,i([ ,i]) = 0. 

(ii) Assume additionally that (QJ) admits an invariant distribution Q. Then for uq G i 1 ([0, 1]), for u the 
corresponding weak solution to (QJ), u(t) goes in law to Q as t — > 00. 
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Clearly, (I) holds if b is C 1 with b' < — e < (choose p(z) = ez) or if a is C 1 with |ct'| > e > (choose 

p(z) — (ez) 2 ). One may also combine conditions on b and a. 

But (X) also holds if 6 is C 1 and if 6' < vanishes reasonably. For example if b(z) = — sg(z) min(|z|, |z| p ) 

for some p > 1, choose p = ep p with e small enough and /0 p (z) = z p for z e [0, 1] and p p (z) = pz — p + 1 

for z > 1. If &(z) = — z — sinz, choose p = 6/33 with e small enough. 

One may also consider the case where a is monotonous with a' vanishing reasonably. 

Let us now compare Theorems [4] and [5] with known results. The works cited below sometimes concern 
different boundary conditions, but we believe this is not important. 

• Sowers [T2] has proved the existence of an invariant distribution supported by C([0, 1]), assuming ("H), 
that a is bounded and that b is of the form b(z) = —az + /(z), for some bounded / and some a > 0. He 
obtained uniqueness of this invariant distribution when er is sufficiently small and bounded from below. 

• Mueller [3] has obtained some surprising coupling results, implying in particular the uniqueness of an 
invariant distribution as well as a the trend to equilibrium. He assumes (H), that a is bounded from 
above and from below and that b is non-increasing, with |6(z) — b(r)\ > a\z — r\ for some a > 0. 

• Cerrai [1] assumed that a is strictly monotonous (it may vanish, but only at one point). 

(i) She obtained an asymptotic confluence result which we do not recall here and concerns, roughly, the 
case b(z) ~ — sg(z)|z| m asz-> ±oo, for some m > 1. 

(ii) Assuming {%), she proved uniqueness of the invariant distribution as well as an asymptotic confluence 
property, under the conditions that for all r < z, b(z) — b(r) < X(z — r), and \a(z) — a(r)\ > p\z — r\, for 
some p > and some A < p 2 /2 (if b is non-increasing, choose A = 0). 

Thus the main advantages of the present paper are that the uniqueness of the invariant measure requires 
very few conditions, and we allow a to vanish (it may be compactly supported). 

Example 1. Assume (H) and that b strictly decreasing. Then there exists at most one invariant dis- 
tribution. If b{z) — — z or b(z) = — z — sinz or b{z) = — sg(z) min(|z|, \z\ p ) for some p > 1, then we 
have asymptotic confluence of solutions. Here to apply [HI [5] one needs to assume additionally that a is 
bounded from above and from below, while to apply [1], one has to suppose that a is strictly monotonous. 

Example 2. Assume {%), that b is non- increasing and that a is strictly monotonous. Then there exists 
at most one invariant distribution. 

If furthermore a is C 1 with < c < a' < C, then we get asymptotic confluence of solutions using [J] 
or Theorem \5\ (here [HI [S] cannot apply, since a vanishes) . But now if a' > reasonably vanishes then 
Theorem \5\ applies, which is not the case of [1]: take e.g. a{z) = sg(z) min(|z|, \z\ p ) for some p > 1, or 
a(z) — z + sinz. 

Example 3. Consider the compactly supported coefficient a(z) = (1 — z 2 )ln z \<i\. Assume that b is C 1 , 
non-increasing, with b'(z) < — e < for z e (— oo, —I) U {0} U (I, +oo). Then Theorems Q] and [5] apply, 
while [HHE] do not. 

Observe here that if b(zo) = for some zq ^ (— 1, 1), then u{t) = zq is the (unique) stationary solution. 
If now b(— 1) > and 6(1) < 0, then the invariant measure Q (that exists due to Sowers [12 ) is unique 
and one may show, using the comparison Theorem of Donati-Pardoux 6 , that Q is supported by [—1, 1]- 
valued continuous functions on [0, 1]. 

However, there are some cases where [121 [4] provide some better results than ours. 

Example J^. If a(z) — pz and b(z) = Az, then u(t) = is an obvious stationary solution. Theorems 2] and 
[5] apply if A < and |A| + \p\ > 0. Cerrai [1] was able to treat the case A > provided p 2 /2 > A. 

Example 5. If a is small enough and bounded from below and if b(z) — —az + h(z), with a > 
and h bounded, then Sowers |12j obtains the uniqueness of the invariant distribution even if b is not 
non-increasing . 
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1.5. Plan of the paper. In the next section, we prove some inequalities concerning the L 1 ([0, l])-norm of 
the difference between any pair of mild solutions to ([!} . Section|3]is dedicated to the proof of our existence 
result Theorem [2j Theorems [4] and [5] are checked in Section [4] We briefly discuss the multi-dimensional 
equation in Section [5] and conclude the paper with an appendix containing technical results. 

2. ON THE i 1 ([0, l])-NORM OF THE DIFFERENCE BETWEEN TWO MILD SOLUTIONS 

All our study is based on the following result. We set sg(z) = 1 for z > and sg(z) = —1 for z < 0. 

Proposition 6. Assume (T-L). For two bounded- measurable initial conditions uq,vq, let u,v be the cor- 
responding mild solutions to ([I]]. Then, enlarging the probability space if necessary, there is a Brownian 
motion (B t )t>o such that a.s., for all t > 0, 

(7) \\u(t)-v(t)\\ L i( [0 ,i]) < IK-Vo||ii([o,i])+ / \W(u(s)) - a(v(s))\\ L mo,x\)dB s 

Jo 

t n i 

sg(u(s, x) — v(s, x))(b(u(s, x)) — b(v(s, x)))dxds. 



o Jo 



Proof. We divide the proof into several steps, following closely the ideas of Donati-Pardoux [6j Theorem 
2.1], to which we refer for technical details. 

Step 1. Consider an orthonormal basis (ek)k>i of L 2 ([0, 1]). For k > 1, we set B^ = J Q J Q ek(x)W(ds, dx). 
Then (B h )k>i is a family of independent Brownian motions. For n > 1, consider the unique adapted 
solution u n eL 2 {9. x [0,T],V), where V = {/ € ^([0, 1]),/' (0) = /'(l) = 0}, to 

u n (t,x) =u (x) + / [d xx u n {s,x)ds + b{u n {s,x))]ds + y2 / a(u n (s, x))e k (x)dB h s . 
Jo k=1 Jo 

We refer to Pardoux [TU] for existence, uniqueness and properties of this solution. We also consider the 
solution v n to the same equation starting from vq. Then, as shown in [B], 

(8) lim sup E[\u n (t,x)-u(t,x)\ 2 + \v n (t,x)-v(t,x)\ 2 } = 0. 

™ [0,T]x[0,l] 

Step 2. For e > 0, we introduce a nonnegative C 2 function e such that </> e (z) = \z\ for \z\ > e, with 
l^eC 2 )! < 1 an d < <f>"(z) < 2e _1 l| z | <e . When applying the Ito formula (see [BJ for details), we get 

(9) / (j),:(u n (t,x) — v n (t,x))dx — / 4> e {uo{x) — vo(x))dx 



Jo 

: r I. 



+ / 4>' e (u n (s,x) — v n {s 1 x))d xx [u n {s,x) — v n (s, x)]dxds 
Jo Jo 

+ 11 (t)' i (u n (s,x)-v n (s,x))[b(u n (s,x))-b{v n (s,x))}dxds 



JO 

t /■! 



+ J2 / <t>'e{u n (s,x)-v n (s,x))[o-(u n (s,x))-o-(v n (s,x))]e k (x)dxdB* 
k=1 Jo Jo 

oE/ / €(u n (s,x)-v n (s,x))[a(u n (s,x))-a(v n (s,x))] 2 e 2 (x)dxd S 



fc=i 



^0 



2 

^.ii+m+im+m+m- 

Since \z\ < (f) e (z) < \z\ + e for all z, wc easily get, a.s., 

lim/ (/> e (w"(£,x) -v n (t,x))dx = \\u n (t) - v n (t)\\ L i , [01]) and lim l\ = ||u - Uo||n([o,i])- 

e->0 Jn e— >0 
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An integration by parts, using that d x [u n (t, 0) — v n (t, 0)] = d x [u n (t, 1) — v n (t, 1)] = shows that 



iHt) 



t rl 



&(u n (s, x) - v n (s, x))[d x (u n (s, x) - v n (s, x))} 2 < 0. 



JO 

2 / rv i 



Since $(z - r)(a(z) - a(r)) 2 < Ce -:l l| z _ r |.< e |z - r\ 2 < Ce by (H), we have 7 e 5 (i) < Cnte, whence 

lim/_ 5 (f) =0 a.s. 

e-K) 

Using that \<f>' e (z) — sg(z)| < l{| z |<e} and (Ji), one obtains a.s. 



lim 



< lim 



7 e(*)-/ / sg{u n {s,x)-v n (s,xj)(b(u n (s,x))~b{v n {s,x)))dxds 
Jo Jo 

t r l 

l|u»(a,a;)-t;»(«,x)|<el&(«"(S)a:)) - &(«™(s, x))|cfa%te < lim Cie = 
o Jo e ~ >0 



Similarly, 



lim E 

£->0 



If W-V / / sg(u n (s,x)-t; n (s,a ; ))[(7(u n ( S ,x))-c7K( s ,a ; ))]e fc (x)da ; dB s A 

fc=1 Jo io 



= 0. 



Thus we can pass to the limit as e — > in © and get, a.s., 
IK(t)-u n (*)IUi([o,i]) < \\ u o ~ vo\\m[oA]) 



(10) 



sg(u n (s, a;) - v n (s, x))[b(u n (s, x)) - b{v n {s, x))]dxds 
Jo 
n „t pi 

+ S / / s g(w"(s,x)-w"(s,a;))[cr(u™(s,a;))-cr(w"(s,a;))]e fe (a;)dxd J B s fc . 



fe=i 



JO 



Step 3. Using (H), there holds, for all r 1 ,z\ 1 r 2 ,z 2 in R, 

(11) |sg(r x -«i)[o-(n) -ct-(«i)] -sg(r 2 - z 2 )[cr(r 2 ) -cr(z 2 )]| < C(|ri -r 2 | + |zi - z 2 \), 

(12) |sg(n - zi)[b(ri) - b(z x )} - sg(r 2 - z 2 )[b(r 2 ) - b(z 2 )]\ < C(\n -r 2 \ + \z x - z 2 \). 

Indeed, it suffices, by symmetry, to check that sg(n — zx)[cr(ri) — &(zi)] — sg(r 2 — zi)[a(r 2 ) — o"(zi)]| < 
C\r\ — r 2 \. If Sg(ri — z\) = sg(r 2 — z 2 ), this is obvious. If now n < Z\ < r 2 (or r\> z\> r 2 ) we get the 
upper-bound |cr(n) + er(r 2 ) - 2ct(zi)| < C(|n - Zi| + \r 2 - Zi\) = C\n - r 2 \. 
Using ([8j) , it is thus routine to make n tend to infinity in (|10[) and to obtain, a.s., 

\\u(t)-v(t)\\ L i {[0A]) < \\u Q - w ||li([o,i]) + / / sg(u(s,x)-v(s,x))[b(u(s,x))-b(v(s,x))]dxds 

Jo Jo 

+ 2 / / sg(u(s,x) — w(s,a;))[c7'(u(s,x)) — ct(v(s, x))]ej;(x)da;<iB s . 



(13) 



For the last term, we used that, by the Plancherel identity, setting for simplicity 



a n (s,x) = sg(w"(s,a;)-i;™(s,a;))[cr(u ,l (s,a;))-cr(w"(s,x))], 
a(s,x) = sg(u(s,x) — v(s,x))[cr(u(s,x)) — a(v(s,x))], 
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there holds 



E 



n r t p l 



oo „t „i 



( V / / a n {s,x)e k {x)dB k s - V / / a{s,x)e k {x)dB k s 



- / E IE(/ {""fe 1 ) -«(s, x)\e k (x)d. 
Jo fc > 1 Jo 

||an(s)-a(s)|||2 ([0il]) ds+ ^ 



ds + E 



fe>n+l 



a(s,x)ek(x)dx 



ds 



k>n+l 



E 



\2i 

a(s,x)ek(x)dx J ds =: /„(£) + J n (i). 



Using (HI]) and then ©, /„(*) < C/ ' J* E[\u n (s,x) - u(s,a;)| 2 + \v n (s,x) - v(s,x)\ 2 }dxds tends to as 
n —$■ oo. Finally, J„(t) tends to because J2k>i Jo "^Klo a ( s , x ) e k(x)dx) 2 ]ds — J Q E[||a(s)||| 2 ([o i])]°^ s < 
C J J E(|u(s, a;) — v(s, x)\ 2 )dxds < oo. 

Step 4- A standard representation argument (see e.g. Revuz-Yor |lll Proposition 3.8 and Theorem 3.9 p 
202-203]) concludes the proof, because the last term on the RHS of (|13[) is a continuous local martingale 
with bracket 

2 



L 2 ([0.1])^ S - 



' E ( / s g( w ( s ' x ) _ w ( s > x )) I '!'"! 5 ) a;)) ~ o-(v(s, x))]efc(x)da; j ds = \\a(u(s))~a(v(s))\\ 

We used here again that (efe)fc>i is an orthonormal basis of L 2 ([0, 1]). □ 

Corollary 7. Adopt the notation and assumptions of Proposition^ For all 7 £ (0, 1), all T > 0, 



E 



7/2' 



sup||«(*)-«(t)||2i ([0|1]) + I y \WHt)) ~ o-(v{t))\\ 2 L2{[01]) dt 



< C*{>,7,t|I m — ^OUi^fo,!])! 



where Cb,-y,T depends only on b,j,T. 

Proof. Let C be the Lipschitz constant of b. Denote by L t the RHS of (O . The Ito formula yields 



\u(t)-v(t)\\ L i {[0A 



,-Ct 



< L,e 



-Ct 



K - w ||li([o,i]) _C< / e 5 i^s 



II 



+ / ||a(u(s))- ( T( W (s))|| L2([ o,i])e- 0s d^ 
Jo 



JO 



-Ca 



sg(u(s, x) — v(s, x))(b(u(s, x)) — b(v(s, x)))dxds. 



But J sg(u(s, a;) — u(s,x))(6(u(s, a;)) — 6(w(s, x)))da; < C||u(s) — v(s)||i,i([o,i]) < CL S . Hence 



\\u{t) - u(t)||ii([o,i])e < ||uq - Uo|Ui([o,i]) + / lk(«M) -0"(«( s ))IU 2 ([o,i]) e s dB s =:M t 



,-Ca, 



E 



Hence M t is a nonnegative local martingale with bracket (M) t = L \\<j(u(s)) — cr(v(s))||| /2( .r ^ 
Applying Lemma [51 we immediately get, for 7 G (0, 1), 

7 /2 

sup ||«(t) - «(t)H2 1([0il]) e- c '" + ( / \\a(u(s)) - v(v(s))\\ 2 L2([01]) e~ 2Cs ds 
[0,oo) \Jo 

< C 7 \\u Q ~v \\l 1([01]) . 
The result easily follows. 
Finally, one can say a little more when b is non-increasing. 



e~ 2Cs ds. 



□ 
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Corollary 8. Adopt the notation and assumptions of Proposition^ and assume that b is non-increasing. 
Then for all 7 € (0,1), 

7 

E 



sup \\u(t) - «(t)||2i ([0)1]) + (J \\b(u(t)) - b(v(t))\\ LH[0A]) dt 



\<r(u(t)) - o-(v(t))\\l Hm) dt 



7/2 
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Li([0,l])' 



where C 7 depends only on 7. 

Proof. Since b is non- increasing, Proposition |6] yields 



u(t)-i>(t)||ii([ ,i])+ / ||6(u(s))-6(w(s))|| L i ([0 ,i])ds 



< ||wo-vo||i,i([o,i])+ / lk(«(s)) -ff(«(s))||i2([o,i])dB 8 =:M t , 
Jo 

which is thus a nonnegative martingale with bracket (M) t = L \\a(u(s)) — o~(v(s))\\ 2 L 2(i ^-.ds. Lemma[5] 
allows us to conclude. □ 



3. Existence theory in i 1 ([0, 1]) 

The goal of this section is to give the 

Proof of Theorem^ We start with point (i). Let thus uq € L 1 ([0, 1]) and consider a sequence of bounded- 
measurable initial conditions (wg)n>i such that \\uq — Uo||z,i([o,il) < 2~". For each n > 1, denote by u n 
the mild solution to |T]) starting from Uq. Using Corollary[7] (with 7 = 1/2), we deduce that a.s., 



E 



ra>l 

which implies that 



1/4' 



sup \\u n+1 (t) - U "(t)||^ 2 ([0il]) + U Mu n+1 (t)) - cr(u n (t))\\l Hm) dt 



< 00, 



E 

n>l 



SUp ||M™ +1 (i) - u"(t)|Ui([0,l]) + lk(w" +1 ) - ^(u")IU 2 ([0,T]x[0,l]) 
[0,T] 



< OO. 



Using some completeness arguments, we deduce that there are some (predictable) processes u and S such 
that a.s., for all T > 0, sup [0T] ||u(i)|| L i([o,i]) + Jo H^WlH^o,:!])* < °° and 

limsup |K*)-u"(*)||i,i([o,i]) = 0, lim||5 , -cr(u™)|| L 2 ([OT]x[oa]) = 0. 

n [ 0]T j • n 

Since ct is Lipschitz-continuous, we deduce from the first equality that lim^ ||ct(u) — c(M n )||^i([o.T]x[o,il) = 
0, while from the second one, lim„ \\S — c( mTI )IIl 1 ([o. t ]x[0- 1 D — 0- Consequently, S = o~(u) a.e. and we 
finally conclude that a.s., 



(14) forallT>0, lim I sup \\u(t) - tt n (t)|| i i ([0> i]) + / ||er(u(t)) - a(u n (t))\\h([o,i])dt 



0. 



It remains to prove that u is a weak solution to (fTJ). We have already seen that u satisfies (J3J). Next, for 
tp e Cf ([0, 1]) with ip'(0) = tp'(l) = 0, for t > 0, we know that a.s., A"'* = B t n ' v for all n > 1, where 

ft /-l 

[u n (s, x)ip"{x) + b(u n (s, x))ip{x)]dxds 



A n,V 



1 ,1 

(p(x)u n (t,x)dx — / yj(x)uQ (x)dx 



Jo 



£«,<P 



JO 



cr(« n (s, a;))tp(x)M / (ds, dx). 
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It directly follows from (|T4l and (H) that a.s., 



lim ^ = 

n— >oo 



tp{x)u(t,x)dx — / ip(x)uo(x)dx 



JO 



[u(s, x)ip"{x) + b(u(s, x))ip(x)]dxds. 



We deduce that Bf := lim„ B"' v exists a.s., and it only remains to check that Bf — Cf a.s., where 
Cf := f L a(u(s, x))ip(x)W(ds, dx) . To this end, consider, for M > 0, the stopping time 

tm = inf { r > ' """' "' ' J - 



l^«WJIlL=([o,i]) da + su P / lk(« n («))lli» 

/0 « Jo 

Using p4p and the dominated convergence Theorem, we see that for each M > 

pt/\T M 

jn,ip (tip |2 



2 ([0il]) eto>M 



limEflH 



tt\TM 



Cr ATM |1=limE 



o 



(a(u(s)) - <7(u n {s))<p\\ L2{[0A]) di 



= 0. 



But we also deduce from (fl4|) that a.s., sup„ L \\a(u n (s) 



|^2(r 01 i)rfs < oo for all T > 0, whence 
lmiM^oo tm = oo a.s. We easily conclude that B^ ,v tends to Cf in probability, whence Bf = Cf 
a.s. 

Point (ii) is easily checked: let (uq)„>i be another sequence of bounded-measurable initial conditions 
converging to uq and let (u n ) n >i be the corresponding sequence of mild solutions to {T]). Then necessarily, 
||uq — 2olU 1 ([o,i]) tends to 0, whence, by Corollary[7l sup[ T ] ||u"(t) — w n (t)||i 1 ([o,i]) tends also to 0, in 
probability. Using (fl4|) . we conclude that supr T i \\u(t) — u n (i)||ii([o,i]) tends to in probability. 

We now prove point (iii). For uq and vq in L 1 ([0, 1]), we consider Uq and Vq bounded-measurable with 
| |uq — uo| |i 1 ([o,i]) + 1 l u o — ^ol |i!([o,i]) < 2~ n . We denote by u,v,u n ,v n the corresponding weak solutions to 
IU). In the proof of (i), we have seen that a.s., lim„ supj 0]T ][||w™(i) — M(t)|| i i([ O ,i]) + ||w n (i)~ l; (0lk 1 ([o,i])] = 
and lim n / [||cr(u n (t)) - cr(u(t))||| 2([0]1]) + ||cr(w n (t)) - <r(v{t))\\ 2 L2 n m]dt = 0. Using the Fatou Lemma 
and Corollary [71 we thus get 

7/2" 
E 



7/2' 



sup||u(t)-t;(t)||2 1([0il]) + \J \Wu(t)) - a(v(t))\\i 2{[0<1]) dt 



< lim inf E 



sup \\u n (t) - v n (t)\\l Hm) + U \\a(u n (t)) - *(v n (t))\\imo,i])dt 



< limmfC 7iT ||Mo - u oll2i([o,i]) = ^y/HK - «o||Ji ([0)1]) - 
Point (iv) is checked similarly. 



□ 



4. Large time behavior 

We now prove the uniqueness of the invariant measure. 

Proof of Theorem^ Consider two invariant distributions Q and Q for ([1]), see Definition [3] Let uq be 
Q-distributed and uq be Q~distributed. Consider the corresponding (stationary) weak solutions u, u to 
(fTj). Applying Theorem [2]- (iv) and the Cauchy-Schwarz inequality, L K s ds < oo a.s., where 

K s := K(u(s),u(s)) = \\b(u(sj) - b(u(s))\\ LHm) + \\a(u(s)) - er(fi(«))||ii {[0>1]) . 

Using Lemma fTUl there is a sequence (t n ) n >i such that K tri tends to in probability. Consider the 
function <j>{r) = r/(l + r) on R+, and define * : L l ([0,l}) X I^flO.l]) i-> M+ as V(f,g) = <f>(K(f,g)). 
Then lim„E[*(u(t„),w(i„))] = lim n E[cj>{K t J] = 0. 

We now apply Lemma [TT1 The space i 1 ([0, 1]) is Polish and for each n > 1, C(u(t n )) — Q and £(u(£„)) = 
Q. The function \1/ is clearly continuous on i 1 ([0, 1]) x L 1 ([0, 1]), (because a, b are Lipschitz-continuous). 
Finally, \&(/, g) > for all f ^ g (because V&(/, ff) = implies that b o f = b o g and a o f = cr o c; a.e., 
whence f = g a.e. since (cr, 6) is injective). Lemma QT] thus yields Q = Q. □ 
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Finally, we give the 

Proof of Theorem{5^ Point (ii) is immediately deduced from point (i). Let thus uq,vo £ L 1 ([0, 1]) be 
fixed and let u,v be the corresponding weak solutions to (U}. We know from (Z), the Jensen inequality 
and Theorem HKiv) that a.s., 

p(\\u(t) - v(t)\\ LH[0A]) )dt < / ||p(|u(t)-t;(t)|)||ii ([ o,i])dt 
Jo 

/>oo 

< J \\\Ku(t))-b(v(t))\ + \a(u(t))-a(v(t))\ 2 \\ Ll{m) dt<oo. 

Using Lemma ITUl one may thus find an increasing sequence (t n )n>i such that p(\\u(t n ) — u(tn)||i 1 ([o,i])) 
tends to in probability, so that ||u(4 n ) — "(irOHLVfo 11) a ^ so tends to in probability (because due to I, 
p is strictly increasing and vanishes only at 0). Next, we use Theorem HKiv) with e.g. 7 = 1/2 to get, 
setting A t = SUp [t]0o) \\u(s) - v(s)\\ L i {[0A]h 



E 



A, 



1/2 



Ft. 



,1/2 



< C\\u(t n ) - w(in)|| L i([ 0il] ) -> in probability. 



We used here that conditionally on j v tn , (u(t n + £,x))t>o,xe[o,i] is a weak solution to ([TJ, starting from 
u(t n ) (with a translated white noise). Thus for any e > 0, using the Markov inequality 



fin 



P [A tn > e] = E [P(A t „ > e\ T t J] < E min (l, e" 1 /^ A^ 2 
which tends to as n — > 00 by dominated convergence. Consequently, as n tends to infinity, 

(15) A tn tends to in probability. 

But a.s. s h- > A s = supr s ^ \\u(t) — i>(t)||.Li([o,i]) is non-increasing, and thus admits a limit as s — > 00, 
which can be only due to (IT51) . D 

5. Toward the multi-dimensional case? 
Consider now a bounded smooth domain Del 1 *, for some d > 2. Consider the (scalar) equation 

(16) d t u(t, x) = Au(t, x) + b(u(t, x)) + a(u(t, x))W(t, x), t>0,xe D, 

with some Neumann boundary condition. Here W(dt,dx) = W(t,x)dtdx is a white noise on [0, 00) x D 
based on dtdx. We assume that a, b : M. m- M. are Lipschitz-continuous. 

It is well known that the mild equation makes no sense in such a case, since even if a(u) is bounded, 
Gt- s (x,y)o-(u(s,y)) does not belong to L 2 ([0,t] x D). The existence of solutions is thus still an open 
problem. See however Walsh [13] when a = 1, b(u) = au and Nualart-Rozovskii [8] when a(u) = u, 
b(u) = au. In these works, the authors manage to define some ad-hoc notion of solutions, using that the 
equations can be solved more or less explicitly. In the literature, one almost always considers the simpler 
case where the noise W is colored, see Da Prato-Zabczyk [5]. 

However the weak form makes sense: a predictable process u = {u(t, x))t>o,xeD is a weak solution if a.s., 

(17) forallT>0, sup \\u(t)\\ L i {D) + f \\<r(u(t))\\% (D) dt < 00 

[0,T] JO 

and if for all function <p E C%(D) (with Neumann conditions on dD), all t > 0, a.s., 

u(t, x)<p(x)dx = / uo(x)<p(x)dx+ / [{u(s,x)Aip(x) + b(u(s,x))}dxds + a(u(s,x))ip(x)W(ds 7 dx)}. 
d Jd Jo Jd 

Assume now that tr(0) = 6(0) = 0. Then v = is a weak solution. Furthermore, the estimate of Theorem 
EKiii) a priori holds. Choosing uq e ^(D) and vq = 0, this would imply (fl7|) . Unfortunately, we are not 
able to make this a priori estimate rigorous. 



STABILITY OF THE STOCHASTIC HEAT EQUATION IN Z, 1 ([0,1]) 



11 



But following the proof of Proposition [5] and Corollary [JJ one can easily check rigorously the following 
result. For (ek)k>i an orthonormal basis of L 2 (D), set B\ — L J D ek(x)W(ds, dx). For u £ L°°{D) and 
n > 1, consider the solution (see Pardoux [ID]) to 

pt " pt 

u n (t,x) = u (x)+ / [d xx u n (s,x) + b(u n (s,x))]ds + ^2 / o-(«™(s,a;))efc(a;)dB*. 



fc=i 



Then if o-(O) = 6(0) = 0, for any 7 G (0, 1), any T > 0, 

„T n 

lL!(_D 



(18) 



[0,T] 



su P |K(t)||2i f ^ + -{ / X)( / ^(«"(*,!c))efc(a:)da:) ds[ 



fc=i 



< C6,7,t||wo|I2i (d) , 



where the constant Cb^.T depends only on 7, T, 6 (the important fact is that it does not depend on n). 
Passing to the limit formally in (fTSj) would yield (fT7j) . Unfortunately, (|18|) is not sufficient to ensure that 
the sequence u n is compact and tends, up to extraction of a subsequence, to a weak solution u to ()16[) . 
But this suggests that, when cr(0) = 6(0) = 0, weak solutions to (|16|) do exist and satisfy (fTTj) . 

6. Appendix 

First, we recall the following results on continuous local martingales. 

Lemma 9. Let (Mt)t>Q be a nonnegative continuous local martingale starting from m £ (0, 00). For all 
7 G (0, 1), there exists a constant C 7 (depending only on 7 J such that 



E 



supM7 + (M)^ 2 



[0,. 



< C^m?. 



Proof. Classically (see e.g. Revuz-Yor [IT] Theorems 1.6 and 1.7 p 181-182] ), enlarging the probability 
space if necessary, there is a standard Brownian motion (3 such that M t = m + 0/m) ■ Denote now by 
r a = inf{£ > 0; /3 t = a}. Since M is nonnegative, we deduce that 

(M)^ < r_ m and sup M t < m + sup j5 s . 

[0,00) [0,T_ m ) 

Thus we just have to prove that E[rZ„J + ^[^m] — C 1 rn i , where S m = supr 0iT _ m ) (3 8 . 

First, for x > 0, P[S m > x] — P[t x < r_ m ] = m/(m + x). As a consequence, since 7 G (0, 1), 



E[S 7 



/ P[Sl > x]dx = / 
Jo Jo 



Ah 



dx = to 7 / Ti—dy = C^m? 

./n l + H 1 ^ 7 



'0 ±-rv 

Next, for t > 0, P[r_ m > i] = Ppnfm t] /3 S > — m]. Recalling that inf[o,t] /3 S has the same law as — v£|/3i| 
we get P[r- m >t] = P[\f3i\ < m/y/i}. Hence 



r V 2 l - 



EfrZ^I = / P[tZ£ > t]dt = / P[|/3i| < ™/i 1/7 ]di = / P[(m/|/9i|)T > t]« = m 7 E [IAI" 7 ] 



D 



This concludes the proof, since E [|/?i|~ 7 ] < 00 for 7 G (0, 1). 

Next, we state a technical result on a.s. converging integrals. 

Lemma 10. Let (K t )t>o be a nonnegative process. Assume that A^ = L K t dt < 00. Then one may 
find a sequence (t n ) n >i increasing to infinity such that K tn tends to in probability as n — > 00. 



Proof. Consider a strictly increasing continuous concave function 
Using the Jensen inequality, we deduce that 



1-4 [0, 1] such that 0(0) = 0. 



1 
TJ 



E[<j)(K s )]ds = E 



1 
TJo 



cj>(K 8 )ds 



< 



TJ, 



K s ds 



< E 



A x 



which tends to as T — > 00 by the dominated convergence Theorem. As a consequence, we may find a 
sequence (t n ) n >\ such that lim n E[0(if{ n )] = 0. The conclusion follows. □ 
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Finally, we prove a technical result on coupling. 

Lemma 11. Consider two probability measures ji,v on a Polish space X. Let \& : X x X i— > R + be 
continuous and assume that ^{x,y) > for all x =/= y. If there is a sequence of X x X-valued random 
variables (X n ,Y n ) n >i such that for all n>l, C(X n ) = /i and L(Y n ) = v and z/lim n E[4 , (X n , Y n )] = 0, 
then \x = v. 

Proof. The sequence of probability measures (C{X n , Y n )) n >i is obviously tight, so up to extraction of a 
subsequence, we may assume that (X n ,Y n ) converges in law, to some (X,Y). Of course, C{X) = /i and 
C(Y) = v. Since \M1 is continuous and bounded, we deduce that E[*(X, F)Al] = lim„ E[^(X n , Y n )Al] = 
0, whence &(X, Y) = a.s. By assumption, this implies that X = Y a.s., so that jj, = v. D 
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